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I. INTRODUCTION 



One of the extensions of electrostatics beyond coulomb interactions is the interaction of 
charged conductors with multiple modes of polarization. The relatively simple problem of 
two charged metal spheres has been studied from a theoretical and pedagogical standpoint 
using image charges [4] [5]. In addition, more complex methods using image charges [6 and 
boundary element methods [7] have been developed for multiple-particle configurations. 

Here we derive a new method to calculate the electrostatic interaction energy of charged 
conducting spheres that utilizes multipole expansions of the potential and charge density. 
This method can be expanded to any number of conductors or to an arbitrary precision 
by including higher order multipoles in the calculation. This approach is motivated by the 
assembly configurations of binary nanoparticle super lattices (BNSLs). Guided by a complex 
combination of entropic and electromagnetic interactions, these lattices can assemble in a 
vast array of different shapes and configurations. In order to better understand the assembly 
of these lattices, a number of experimental studies have tried to isolate and analyze the 
different interactions driving the assembly [8j. Here we focus on analyzing the electrostatic 
interactions that guide the assembly of the metal nanoparticles at the center of these lattices. 

Normally in electrostatic problems, most of the higher order interactions can be ignored 
since the monopole and dipole terms dominate the energy equation for large separations. 
Once the scale of the configuration approaches the scale of the nanoparticles, the quadrupole 
or octupole interactions can be on the same order as some monopole interactions [9]. Since 
these assemblies can also include nanoparticles in contact, there is no assumption that can 
justify ignoring multipole terms beyond the dipole. This multipole expansion method allows 
the interaction between any number of conducting spheres to be calculated using an arbitrary 
number of multipoles. This will provide insight into why BNSLs only assemble with certain 
configurations of metal nanoparticles by allowing a direct energy comparison across different 
configurations. Moreover, it will allow us to analyze how the geometry of a configuration 
influences the relative strength of various multipoles in a configuration. 

We proceed in deriving the equations for the energy of the configuration by first looking 
at the energy between two charged spheres. To make this calculation easier, we begin 
by assuming that the two spheres lie on the z-axis of our arbitrarily selected coordinate 
system. Once we find the energy equations in this simplified aligned case, we derive the 

2 



transformation to express the energy equations in any arbitrary frame. The accuracy of this 
method is then checked against existing literature calculations, and then we generalize the 
energy equations to include an arbitrary number of spheres. Finally, we present calculations 
for a number configurations and comment on the accuracy and insights of the multipole 
expansion method. 



II. METHOD AND DERIVATION 



To calculate the multipole expansion of the energy, we start with the energy equation: 

E=- [ pVd 3 x (1) 
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This can be broken up into the pairwise sphere interactions for N spheres: 

N N 



E 



^EE / ( 2 ) 



zl ' =1 Sphi 



where pi is the charge density on sphere i, Vj is the potential on sphere j, and the integration 
is carried over the volume of sphere i. The terms where i — j represent the energy cost of 
holding a charge distribution on the sphere, further referred to as the self-energy. The other 
terms represent the interaction energy, or the energy associated with charge distributions 
and potentials on different spheres. We now find the multipole expansion of the charge 
distribution and potential for a general sphere S. The charge distribution on sphere S (with 
radius r$) is given by: 

Ps = J2srYe m *(0,mr-rs) (3) 

£,m 

where S™*, Y™*{Q, 0), etc. represent the complex conjugates of SJ 1 , Y™(9, 0), etc. The S™'s 
are related to the multipole moments, % m 's, as defined by Jackson (Gaussian or cgs units) 
fTUl p. 145] where: 

Qim 1 



r l+2 r £+2 



sr = ^ = ^r- 2 Yr(e, typsWd'x (4) 



To calculate the potential, we expand the formula for the potential in terms of the sphere's 
multipoles: 
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FIG. 1. Diagram of the two spheres in the aligned case. 
A. Z-Axis Aligned Case 

To calculate the energy explicitly, we start by just considering two spheres, A and B, 
with their charge distributions and potentials given by replacing the label S with labels A 
and B in equations [3] and [5] . Thus the total interaction energy is given by: 



E 



int 



p A V B d 3 x+^ j p B V A d 6 x 



E AB + E BA 



(6) 
(7) 



We will continue by deriving the formula for E AB . The other half of the interaction energy, 
E BA , can be derived by using p B and Va in place of pa and Vb[I]. For the calculation of 
two aligned spheres, the radius term in the potential (7Z in eq. [5]) is the distance between a 
point on sphere A and the center of sphere B. This gives the E AB interaction energy as: 
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8{r - r A )d 3 x (8) 



We can find the form for 1Z using the law of cosines: 



U = V(Rab) 2 + (r A ) 2 - 2(R AB )(r A ) cos# 



+1 



(9) 



Plugging this in and integrating over the radius variable, gives the energy as: 

- 2 ~ /+2 r Y i m *(9,(f ) )Y.T(0',4>') 



f ab V 2tt m >r\r^ f YrWYfV,*) ( . 
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The expression contained inside the radical is the generating function for the Gegenbauer 
Polynomials [H] 

00 

Vl - 2XZ + Z 2 = ^ C n ] ( X ) ZU , \ Z \< 1 (11) 
n=0 

with 2 = , x = cos#, and ( = ^j-k This allows us to express the denominator as a 
power series in (va/Rab) with the coefficients being the Gegenbauer polynomials(Cn^(x)). 
Through a transformation matrix of coefficients, the denominator is changed from an ex- 
pansion of Gegenbauer polynomials into an expansion of spherical harmonics. This is done 
through the following transformation: 

n n 

C<0(cos<9) = ]T4cos fc (£) = J2 c C nk h kp P p (cos0) (12) 

k=0 k,p=0 

= E4v/^^M) (13) 

n 

= E 7 iW*) ^ 

where we defined 7^ as: 



, 2p+l 
k=o T ^ 



U = E4W^rr as) 



and utilized the following identity: 



n 



h kp P p (cos 6) = (cos6) k (16) 

p=0 

where P p (cos9) is the p-th Legendre polynomial [HJ P- 798]. 

The formula for 9' is derived from the law of sines: sm - = (see figure [l]). As 
before, we can expand 1Z using the law of cosines and the Gegenbauer series expansion: 

d ■ I r A sin# \ 

U = 7T — arcsm ^^^^^^^=^^^^^^^= (17) 

l^V 1+ S )2 - 2 S )cos v 

(ST = tt - aresin E^ ( 18 ) 



5 



(I) 

where C t 2 is the t-th Gegenbauer polynomial. Because of the azimuthal symmetry of the 
set-up we see that (j) f = (j). Plugging this all into the energy equation (leaving 9 f in its place 
for simplicity): 

oo n 



,m,m' 72=0 P=° 

CO Up. 

e ab = E Q w 'Dr )n Ii^ ^^^'^'^^ (20) 
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where A^'/ m is defined as: 

t^m,m 'A'B Am* om /oi \ 

*V = 2? + 1 ^ B ' (21) 

^AB 

and C = ^j-k We can also simplify the integral to be performed by utilizing the Wigner 3-j 
symbols to decompose the product of spherical harmonics into a sum. This formula is given 
by [12, p.1057]: 

ir E V (2£+1)(2p 4 +1)(2L+1) f^ L )™( tpL ) m 
= E wfrYrm (23) 

L=\e- P \ 

( £ P l\ 

where I J is the appropriate Wigner 3-j symbol with the convention used in [T21 

\£ z Pz L Z J 

p. 1053] and f2 = (6, (f>). This gives us a final expression for E AB : 

oo n i+p „ 

e ab = £ ^/^(^r^Tj £ w i? Yr(0,<t>)Y?\e\<t>)dn (24) 

*/,m,m' ' n=0 AB P=° L=\£-p\ J 

This is the interaction energy that we will use to calculate the energy for a configuration 
of spheres. To give an example of what the individual energy terms are, we integrate this 
formula explicitly for the first few multipoles. For simplicity in reporting the expansion 
below, we assume that the spheres are of equal size (ta — tb — t). In addition, we will 
assume the two spheres are separated by a distance much larger than their individual radii 
(#A£ << l)- This allows us to take the following simplification [2]: 

ta sin 6 



9 = it — arcsin 



Rab 
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(a)Three dimensional perspective (b)x-z plane 

FIG. 2. Non-aligned configuration with Euler angles a and (3 marked. 



Using this simplification we expand E AB and report the terms of order r 3 (-j^) 3 or lower: 

eAB _ 2 7T r 4 A° B$ 2nr 5 A\B% _2nr 5 A° B° _ 4 tt r 6 

Rab + V3(i? AB ) 2 " V3(R AB y ~ 3(Rab) 3 
2nr 4 A° 2 B° 2 7rrM°i? 2 

^(RabY V$(Rab) 3 
These are the terms of the multipole expansion of the partial interaction energy, E AB , up to 
the dipole-dipole interaction in the special case where the spheres are far apart and aligned 
along the 9 = axis. 



B. Non- Aligned Case 

To generalize for spheres not aligned along the 9 = axis, the expression for 1Z becomes: 

= V(Rab) 2 + (r A ) 2 ~ 2(R AB )(r A ) cosry"^ 

oo n 

= Ra ( b +1) E(^) n ErcStffo r)] (27) 

where 77 and r are the polar and azimuthal angles between and the vector 1Z. When 
the spheres were aligned along the z = axis, we were able to simplify the equations since 
rj = 9 and r = (f). 

We start with same Gegenbauer expansion method presented above, but we also need to 
rotate the spherical harmonics in (77, r) to spherical harmonics in (0, (j)) by using the Euler 
angles (a, /5, 7) that separate the two frames. In our specific case, we are only rotating 
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the center point of a sphere to the z-axis without changing the orientation. This means 
that one of the Euler angles, 7, will not have any effect on the orientation of the spheres. 
Thus we can take 7 = and just use two Euler angles to describe the rotation. The 
following transformation converts the spherical harmonics with (77, r) arguments to spherical 
harmonics with (0, (j)) arguments: 

£ 

Y f m (v, r) = J2 P M, m (0, -P, -a, ) Y e M (6, 0) (28) 

M=-£ 

where T> £ Mm is the standard finite rotations operator that can be found in Landau and 
Lifshitz [Bl P-218]. 

This transformation allows us to easily calculate the interaction energy equation for a 
configuration of multiple spheres. In a configuration with more than two spheres, each 



pairwise interaction energy has the same form as equation [24] when the coordinate systems 
are aligned along each line of action (i.e. along Rab as 77 is in Figure 2b.). Instead of 
recalculating the entire energy series for every orientation, we can utilize the fact that 



the rotation of the spherical harmonics is a linear transformation (see Appendix A). The 
transformation allows us to take each multipole referred in these local pairwise frames and 
transform them all to the same configuration-wide (0, (j)) coordinate system. This expands 
the A™ and B™' multipoles in (77, r) in terms of new multipoles A™ and B™' in the (0,0) 
frame. This gives us the transformation: 

=^Af=fl V mJ°> ~P> ) A f (29) 

M=-£ 

The rotation coefficients can be further decomposed into their three constituent parts: 

^, m (0, ~P, -a, ) = e iM -°d* M>m (-f3)e- iM * (30) 

A? = E SmJtPW*-™* (31) 

M=-£ 

Af* = E d^ m {-P)Af*e iM « (32) 

M=-£ 

Applying this transformation to the multipoles changes the expression for /C to: 
2tt r 2 A 2 / A . \ I * 

r m,m z/r ' A B | i£ ( o\ jM* JMa \ | i£ ( o\ nM' -iM 

~ Of' + 1 p^+l 2^ rf M,m(-Pj^ 6 ^ d M>,m\-P) B i> e 



(33) 



This new form for JC is then plugged into equation [24} This gives us a completely general 
expression for the interaction energy between two arbitrarily oriented spheres. To see what 
the first few energy terms are, we integrate equation [24] explicitly and obtain the formula 
for the interaction energy between the two spheres. Due to limited space we show only the 
monopole-monopole and dipole-monopole terms in the -^-^ << 1 limit below: 



2tt r 2 r 2 A B° 2 \h cos (!) sin (^) r A r%A^ u B^e~ ia 

j^AB _ Z7T r A r B /l n + V 3 V2/ U 

Rab (Rab) 2 



2 7r(cos(f) 2 -sin(f) 2 )ri^?^o° 2^|tt cos (f ) sin (f r|^*5 V 



V3(R AB ) 2 (Rab) 2 



C. Self-Energy 



The self-energy formulas are much easier to derive. Starting with the expression for the 
energy in equation [2j we set i — j and use the charge density and potential definitions from 
equations [3] and [5| This gives the self-energy for sphere A is: 



(35) 



£,£ ,m,m' 



Simplifying and utilizing the orthogonality of the spherical harmonics to perform the inte- 
gration yields: 

ESdf = E ^-/aATAT (36) 

Applying the same process to sphere B yields the same exact result in terms of the £>^'s 
instead of the A™'s. By combining the interaction energies and the self energy from both 
spheres we obtain a energy formula that can be expanded to any order interaction. 



D. Green's Theorem Simplification 

Due to the many summations in the interaction energy formula, the computation of 
relatively small order interactions (like the monopole-dipole interaction) can be extremely 
complicated and asymmetric. For example, the computation time needed to calculate E^_ d 
(a monopole on sphere A and a dipole on sphere B)[3j is much longer than the computation 
time needed to calculate E^ m even though both represent the same electrostatic interaction. 
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Looking at eq. 10, it is clear that the t! exponent makes E AB more difficult to calculate as 
the multipole order increases. 

This calculation can be simplified by making use of Green's reciprocation theorem to 



reduce the order of the exponent in eq. 10 fTUl p. 52]. When constructing the interaction 
energy E AB , we can split the calculation up into each of the different components: 

E AB = Ei B _ m + El B m + Ei B _ d + Ej* d + ... (37) 

From Green's reciprocation theorem, we know that: 

rpAB T^BA (no\ 

In order to calculate E AB _ d , we would need to expand Y™' {0* , (j)) in another series. Thus, it 
is advantageous to use E BA m instead since it avoids having to use any further expansions or 
unnecessary computation time. This renders the interaction energy as: 

E AB = Ei B _ m + Et B m + E B _ A m + + ... (39) 

This same idea is then applied to further terms in the expansion (e.g. replacing Ef B q with 
Ef- A d ). 

In addition, Green's reciprocal theorem also states that E AB = E BA . With proper atten- 
tion to the conjugation of multipoles, there is no need to calculate E BA . Since substitutions 



can rotate our aligned energy into any rotated frame (see section II B), no calculations be- 
yond E AB in the aligned case are needed to obtain the full interaction energy (E int ), for any 
orientation. 



E. Checking the formula 

Here we will check specific terms of the energy equation against other known equations. 
We will use the simpler aligned version of the energy equation for these calculations. The 
first and easiest check is that the monopole-monopole interaction gives the correct coulomb 
energy. For two charged spheres, we use the definitions of the multipole moments in Jackson 
m P- 146], A° = ^= to find: 

F int =F AB _|_ rpBA = ^AB = 4^rjr|<gj = 4 71 T A j (ggg) ( gggj = g A q B 

m m—m * m—m m—m jd td r> 

EtAB EtAB EtAB 

(40) 
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which agrees with the coulomb energy of two point charges. The dipole-monopole and the 
dipole-dipole energies derived from our equations also agree with the Jackson formulas. To 



check the dipole-monopole energy, we use the simplification of section |IID| and plug in the 
values of the multipoles to find: 

rnint T?AB , t^BA ^t^AB ( A1\ 

4 7rr 3 r 2/i0D0 4 7T r\ r\{— ) ( 2 ) 

~ Mrf " vf(w - ^ <w - p - 4B,isM (42) 

where is the dipole on sphere A, and Efidd is the electric field due to a monopole on 
sphere B. This agrees with the formula for the energy of a perfect dipole in an electric field. 
For the dipole-dipole interaction, we start with the formula for the interaction between two 
dipoles given by Jackson [10, p. 147]: 

j? _ PrP2-3(n-pi)(n-p 2 ) fA ^ 

Since this formula is for the energy between two perfect dipoles, we recalculate the dipole- 
dipole term of our interaction energy assuming << 1 in order to compare it to the 



Jackson formula. Using the Green's theorem simplification from section |IID| again and 
taking the dipole-dipole term from the E AB expansion yields: 

jpint t?AB i j?BA 0T? AB 8 7ir\r%A^B^ 1 5 

h d-d = h d-d + h d-d = Zh d-d = — wtb — \3 — + °yy~B — ) ) I 44 ) 

oyriAB) K AB 

Now looking at the Jackson formula, we note that the distance between the two dipoles 
will be Rab- Since the dipoles could either be aligned in the same direction or in opposing 
directions, there is a plus or minus factor in front of the formula: 



mt ±(p lP2 - 3p lP2 ) ±2p lP2 ±2 (y 1 A -ir\A\){^-fr%Bl) ±8^3^030 

d ~ d R 3 R 3 (Rab) 3 HRab) 3 [ ' 

which agrees with our lowest order term from the dipole-dipole interaction energy expansion. 
We can also compare a numerical result from our calculations with an analytical solution 

done by Maxwell[14, p. 275]. Maxwell finds that two touching spheres of radius r held at 

potential V = 1 statvolt each have a charge of Q = rlog(2). Taking r = 1 cm and using the 

equations for the energy stored in a capacitor we find: 

c= Q^a L = 21o|(2) =2lQg(2) (46) 

11 



TABLE I. Comparison of Calculated Energies with Analytic Maxwell Energies 



Highest Order Term Energy (ergs) Percent Difference 
Dipole-Dipole 1.450672 0.553 % 

Quadrupole-Quadrupole 1.444552 0.128 % 

Octupole-Octupole 1.443789 0.076 % 



= ^ = i^T '- 442695 ergs < 47) 

where E cap is the energy stored in the capacitor. This is extremely close to our energy 
values for the Rab = 2 limiting case. In addition, our formula converges rather quickly 
to the correct value. To calculate the energies reported above, we explicitly performed 
the integrations from our final energy equation using Mathematica. We then minimized the 
energy equation utilizing the matrix formalism detailed in the next section. The Gegenbauer 
expansion of the denominator, 1Z ) was taken out to the sixth Gegenbauer polynomial. This 
was the point at which we started getting marginal returns in the percent difference between 
our calculated energy and Maxwell's figures. The second Gegenbauer expansion used in the 
expansion of 0* was taken out to the first term. This provided improvement over the small 
angle approximation and left the integrals still analytically computable by Mathematica. 



F. Generalizing for more spheres 

These formulas are generalized for any number of spheres by writing the multipoles as a 
vector and the individual terms of the energy as the entries of a matrix. This gives us the 
following matrix equation for the energy between any number of spheres: 

E = ^M+KM , with M = ^Aq A\ . . . B° . . . C ° . . .) T (48) 

where M is a vector containing all the multipoles for all the spheres to a particular order 
and K is the energy matrix that contains all the coefficients for the different multipole 
interactions. The diagonal terms of K are the self-energy terms, while generally, K n5m is the 
interaction between the M n and the M m multipoles. The left section of K for a three sphere 
calculation up to dipoles is produced in the appendix. The minimum of equation [48] occurs 
in the trivial case when all the multipoles are equal to zero. To examine the non-trivial cases 
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we need to find the constrained minimum by minimizing the function: 

T = ^M»KM - ^ Wi ■ M - Qi) 



(49) 



where are the undetermined multipliers for the constraints. The term q is a dummy 
variable that uses the dot product q • M to select out a specific multipole to constrain. The 
variable Qi is the constrained value for that specific multipole. For example, to constrain 
the monopole term on sphere A to be one, the second term of the equation would take the 
following form: 



A(c-M-Q) = [ 







V : J 



KK - 1) 



(50) 



While we can arbitrarily constrain any multipole, our calculations only constrain the 
monopole terms since they correspond to the total charge on each sphere. For a configuration 
with three spheres, there are three different constraint equations; one for the monopole term 
on each sphere. 



In order to minimize J 7 , we take the derivative of equation 49 and set it equal to zero. The 
derivative is easy to calculate since K is a Hermitian matrix, and this yields the multipoles, 
M, that minimize the constrained equation: 

V T = KM - \ci = 

i 

=4> M = K" 1 ^A i c i 

i 

By using this matrix formalism and solving for the unconstrained multipliers we can 
find the multipoles that minimize the electrostatic energy for any given number of charged 
conducting spheres. 



III. DATA - ENERGY OF SPECIFIC CONFIGURATIONS 

Here we present the configuration energy and the work required to remove a bounded 
sphere for various configurations. The configuration energy includes both the interaction 
and self energies of all the spheres except the monopole self energy. Since we constrain the 
monopole terms to constant values, we omit them from the energy data. The work required 
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to remove a sphere is calculated by finding the difference in configuration energy between 
the original configuration and the configuration with one sphere removed. In the case of two 
spheres, the work needed to remove a sphere is the same as the configuration energy. In the 
case of three spheres, we subtract the total configuration energy from the energy of the two- 
sphere system that results after we remove a sphere. For asymmetrical configurations, we 
compare the work required to remove each of the different spheres and report the minimum. 

The spheres in the figures are shaded with their charge density. The coloring goes from 
blue to red as the charge density goes from negative to positive (color online, dark to light 
in grayscale). In addition there are lines drawn around the z-axis of each sphere to help 
qualitatively show the orientation of the charge on the spheres. The figures are meant to 
show alignment of the various higher-order multipoles relative to the other spheres in its own 
figure. The shading scale was varied between different figures to provide enough contrast 
within a figure. 

In constructing these configurations, we label spheres with letters (A,B,C.) usually 
placing sphere A at the origin as a convention. Since the equations are rotationally invariant 
usually only one specific configuration is present (i.e. no more than one equilateral triangle 
with one charged sphere will be looked at since all three possibilities are equivalent). 

For the calculations, we used the same parameters discussed in section [HE] when we com- 
pared our energy values with Maxwell's calculations. Charged spheres had their monopole 
constrained so that the sphere has one statcoulomb of charge. All spheres have radius of 1 
cm, and all energies are reported in ergs (Gaussian or cgs units). 



A. Two Sphere Configurations 

1. Discussion of Two Sphere Configurations 

In the case of two spheres, there are only two different configurations of interest. The 
first is when both spheres have equal charges. This induces multipoles equal in magnitude 
for both spheres and equal binding energies for each sphere. When only one sphere is 
charged, the higher order multipole terms have a much stronger role. Looking back at the 
explicit interaction energy expansion, the first three terms of the energy expansion are the 
monopole-monopole energy and the two monopole-dipole interactions. However, only one of 
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TABLE II. Two spheres Q 



(a)Both Charged 



Highest Order Interaction 


Energy (ergs) 


Octupole-Octupole 




0.443789 


Quadrupole-Quadrupole 




0.444552 


Dipole-Dipole 




0.450672 


(b)One charged 


Highest Order Interaction 


Energy (ergs) 


Octupole-Octupole 




-0.0557567 


Quadrupole-Quadrupole 




-0.048093 


Dipole-Dipole 




-0.0336496 



The energies differ by 1 compared with the energies cited in the Maxwell's calculation since we do not 
include the self-energy due to the constrained monopole. 



s 



: 



(a) Both charged (b) Bottom sphere charged 

FIG. 3. Comparison of the two configurations with the total charge distribution shaded on the 
surface. 

those three terms is non-zero when one sphere is left uncharged. This results in the binding 
energy for the more weakly bound sphere (Sphere B) to increase 43% when the expansion 
is extended from the dipole-dipole interaction up to the quadrupole-quadrupole interaction. 
If we take the calculation up to the octupole-octupole interaction, the binding energy of 
that sphere increases by another 15%. This is compared to only about a 0.2% difference 
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(a) Both charged (b) Bottom sphere charged 

FIG. 4. Comparison of the two configurations with only the quadrupole charge density shaded on 
the surface. 



between the quadrupole-quadrupole and octupole-octupole calculations in the configuration 
with both spheres charged. If we had arbitrarily cut off the calculation at the dipole-dipole 
term of the configuration with only one sphere charged, we would have seen a binding energy 
for sphere B that is about one-third of the strength it actually is. This shows evidence for 
why higher ordered terms cannot be arbitrarily ignored in these energy calculations. Also, 
since the configuration energy relies heavily on higher order multipoles, the energy values 
do not converge quickly when only one sphere is charged. 
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TABLE III. Three spheres 
(a)Total configuration energies vs. Highest order interaction (energy in ergs) 





A 


B 


C 


D 


E 


Octupole-Octupole 


-0.111602 


-0.0999026 


-0.09723 


-0.0702905 


-0.0690104 


Quadrupole-Quadrupole 


-0.0965803 


-0.0871968 


-0.0889388 


-0.0570036 


-0.0602305 


Dipole-Dipole 


-0.0596624 


-0.0627721 


-0.0605811 


-0.0438303 


-0.0400716 


(b)Work to remove weakest bound sphere 


to infinity vs. 


Highest order interaction (enerj 


iy in ergs) 




A 


B 


C 


D 


E 


Octupole-Octupole 


0.055845 


0.0441459 


0.0414733 


0.0145338 


0.0132537 


Quadrupole-Quadrupole 


0.0484873 


0.039103 


0.0408458 


0.00891062 


0.0121374 


Dipole-Dipole 


0.0260128 


0.0291225 


0.0269315 


0.0101807 


0.00642202 



B. Three Sphere Case 



u 



(a) Equilateral triangle (config. A) (b) Right triangle, charged apex 

(config. B) 

FIG. 6. Comparison of two configurations with the total charge configuration shaded on the surface 
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(a) Equilateral triangle (config. A) (b)Right triangle, charged apex 

(config. B) 



FIG. 7. Comparison of the two configurations with only the quadrupole charge densities shaded 
on the surface 



For the three sphere case, the minimal energy configuration is the equilateral triangle 
(configuration A). This is followed by the isosceles triangle (config. B) and the middle- 
charged line of spheres (config. C). As more multipoles are added, the minimal energy 
configuration changes. When only dipoles are taken into account, both the isosceles triangle 
and the middle-charged line of spheres are favored over the equilateral triangle. This is due 
to minimal dipole-dipole interaction between the two uncharged spheres in each configura- 
tion. The isosceles triangle configuration induces mutually orthogonal dipoles on the two 
uncharged spheres, reducing the energetic cost of the configuration. In the middle-charged 
line of spheres, the monopole on the middle sphere dominates the energy equation through 
the monopole-dipole terms. The antiparallel dipoles that are induced on the two uncharged 
spheres have only a slight effect on the configuration energy. However, this effect is strong 
enough to give the configuration a slightly higher energy than the isosceles triangle. 

Once we add quadrupoles, the dipole advantages discussed above are much less influen- 
tial, causing the equilateral triangle becomes the minimal energy configuration. Since all 
the dipoles in the equilateral triangle line up to the same point in space, there is much more 
repulsive interaction between the dipoles than in the isosceles triangle configuration. This 
caused it to have a misleadingly higher configuration energy when only dipoles are consid- 
ered. With the quadrupoles added, the equilateral triangle experiences the largest increase 
in binding energy, making it the lowest energy configuration. The symmetry of this con- 
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figuration causes the energy to be more dependent on quadrupole terms than the isosceles 
triangle or the middle-charged line of spheres. This shows that some configurations show 
non-regular convergence when different multipoles are added to the calculation. In particu- 
lar, the work to remove a sphere from configuration D is not even monotonically decreasing 
as we add higher order terms, demonstrating the significant effect that the geometry of the 
configuration has on the binding energies. 



In addition, configurations with uncharged spheres equidistant from the charged sphere 
have stronger binding energies. The line with a charge on the end (config. D) and the 
isosceles triangle with a charge on one of the legs (config E) are the most weakly bound 
configurations. Both of these configurations are asymmetrical in the sense that the two 
uncharged spheres are not equidistant from the charged sphere. Judging from the charge 
distribution diagrams of these configurations, it seems that the asymmetric configurations 
restrict the polarization of one of the uncharged spheres. This prevents the configuration 
energy from lowering as far as the other three configurations that all feature the uncharged 
spheres equidistant from the charged sphere. 



We also see that the energy convergence of the equilateral triangle's configuration energy is 
similar to the the two sphere case. The binding energy increases by 86% when the calculation 
order is taken from dipoles up to quadrupoles and a further 16% when octupole interactions 
are added in, comparable to the convergence seen in the two sphere case. However, there 
is variability in energy convergence for different configurations. For example, the binding 
energy for configurations C and D only increase by 1.5% and 8% respectively when the 
calculation is extended from quadrupoles to octupoles. 
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TABLE IV. Four spheres 
(a)Total configuration energies vs. Highest order interaction (energy in ergs) 





A 


B 


C 


D 


E 


Octupole-Octupole 


-0.169347 


-0.156567 


-0.135689 


-0.129386 


-0.112597 


Quadrupole-Quadrupole 


-0.148192 


-0.13541 


-0.119614 


-0.11312 


-0.099499 


Dipole-Dipole 


-0.0804022 


-0.0822249 


-0.0882232 


-0.068025 


-0.0672607 


(b)Work to remove weakest bound sphere 


to infinity vs. 


Highest order interaction (enerj 


iy in ergs) 




A 


B 


C 


D 


E 


Octupole-Octupole 


0.0577453 


0.0449653 


0.0357864 


0.0177843 


0.0126944 


Quadrupole-Quadrupole 


0.0516117 


0.0388297 


0.0324172 


0.0165397 


0.0123022 


Dipole-Dipole 


0.0207398 


0.0225625 


0.0254511 


0.00836263 


0.00448864 



C. Four Sphere Configurations 




(a) Blunt end charged (config B) (b) Sharp end charged (config D) 

FIG. 9. Comparison of two rhombus configurations with the total charge configuration shaded on 
the surface 
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(a) Blunt end charged (config B) (b) Sharp end charged (config D) 

FIG. 10. Comparison of the two configurations with only the quadrupole charge densities shaded 
on the surface 

For configurations with four spheres, the minimal energy configuration is the regular 
tetrahedron (config. A). Much like in the three-sphere case, the minimal energy configuration 
changes as additional multipole terms are added to the calculation. When only dipoles 
are considered, the minimal energy configuration is the non-regular tetrahedron (config. 
C). Since most of the faces of this tetrahedron are isosceles triangles, the behavior of this 
configuration is similar to the isosceles triangle that we analyzed in the three-sphere case. As 
before, the configurations that induce dipoles in mutually orthogonal direction will appear 
to have an energetic advantage when the calculation only includes dipole terms. 

The regular tetrahedron (config. A) mimics the behavior of the equilateral triangle, 
with all of the spheres' induced dipoles focusing on a common point. This causes strong 
interactions between the induced dipoles and makes configuration A look less favorable when 
only including dipole interactions in the calculation. Thus, when only dipoles are included in 
the calculations, there is an intrinsic bias for configurations that feature mutually orthogonal 
dipoles. This bias, caused by arbitrarily cutting off the calculation at the dipole-dipole 
interaction, results in incorrect conclusions for the minimal energy configuration. 

Once quadrupoles are added, the configurations take the same order as when the cal- 
culation is carried out to octupoles. The regular tetrahedron (config. A) is the preferred 
configuration with the charged blunt-end rhombus as the next most favored. Note the signif- 
icant energy difference between the two rhombus configurations (configs. B and D). When 
the blunt-end particle is charged in a rhombus, the charged sphere is in contact with all of 
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the other three spheres in the configuration causing significant polarizations in the three un- 
charged sphere. However, in the sharp-end charged rhombus (config. D), the charged sphere 
is only in contact with two other spheres. This causes a much smaller polarizing effect on 
the third sphere and makes the configuration less energetically favorable when compared 
with the blunt-end charged rhombus (config B). 

IV. CONCLUSION 

Here we presented a new method utilizing multipole expansions to calculate the energy 
of charged conducting sphere configurations. This method can be used for any number of 
spheres in an arbitrary configuration, and it is easily expanded to any multipole order. 

As seen from the configurations we considered, cutting off the energy expansion at ar- 
bitrary multipole orders (like at the dipole-dipole or quadrupole-quadrupole level) gives an 
incorrect picture of the energetically preferred configurations. In addition, there are intrinsic 
biases towards particular configurations. Configurations that induce mutually orthogonal 
dipoles, were favored when we cut off the energy expansion at the dipole-dipole interaction. 

In this study, we truncated the energy expansion at the octupole-octupole level. In 
the case of two charged spheres, this was sufficient since the energy converged quickly; the 
energy only increased 0.2% when the calculation was taken from quadrupoles up to octupoles. 
However, this is not the case when we move on to more general configurations with multiple 
spheres. In general, the more spheres in a configuration, the more the energy calculations 
rely on higher order multipoles, and the less quickly the energy converges. The calculations 
here were also subject to an arbitrary cut-off for interactions at the octupole-octupole level. 
Given more computing power or patience, the calculations could be extended even further. 

We can extend this study by analyzing further how different multipoles interact in various 
configurations. For example, a family of configurations could be found where the configu- 
ration energy only slightly changes with the addition of quadrupoles. This implies that the 
specific geometry of the configuration reduces or neutralizes the quadrupole-quadrupole in- 
teractions. By identifying similar configurations, more complex tailor-made materials can be 
created that utilize configuration geometry to reduce individual inter-particle interactions. 

This paper was the product of an independent research project done by Alex Moore under 
the supervision of Thomas Witten while at the University of Chicago. In addition Nathan 
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Krapf, also at the University of Chicago, helped to extensively revise this paper. Alex would 
like to thank both Prof. Witten and Nathan Krapf for their patience and guidance while 
working on the project, as well as Dmitri Talapin for his support and for providing the initial 
motivation for this study. 



Appendix A: Rotating the Multipoles is Equivalent to Rotating the Spherical Har- 
monics 

Proof. The charge distribution on a sphere and the energy between two spheres are invariant 
with respect to any choice of coordinates. Therefore, given two systems of coordinates, 
X = {#,0} and x = a multipole in the x frame can be expanded in terms of the 

multipoles in the x frame for any point x = 0}, (/)}: 

t 

A?Y e m (e,<t>)= ^ M (M) (Al) 

M=-£ 

Now we can rotate the Y £ M (9, 0) into the x frame by utilizing the identity [13, p. 218]: 

Y e M (9j) = £ [v^Mia^^Y^eA)] (A2) 

N=—£ 

where (a, /3, 7) are the Euler angles needed to rotate the x coordinate system to the x 
coordinate system . Plugging this in: 

£ t 

AfYnOA) = £ Af ( J2 ^,m(«,A7)^(^)) (A3) 

M=-£ N=-£ 
£ £ 

at Y e m (e, <f>)=J2 Y t N ( > M E v ku{^ P, 7)) (A4) 

N=—£ M=-£ 

By the linear independence of the Y™ 's, we must have N = m. This removes the summation 
over N and presents us with our final indent ity: 

£ 

=^AT= J2 2>m,MKA7)Af (A5) 

M=-£ 

If we are trying to rotate a conjugated multipole, we just use the rules of complex conjugation 
of sums and products to get the identity: 

£ 

^A?*=Y1 (^,m(«> A 7))*Af * (A6) 
M=-e 

□ 
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TABLE V. Configuration Details. Location of charged sphere in bold, 
(a) Three sphere configurations 



Configuration 


Letter Sphere A 


Sphere B 


Sphere C 


Sphere D 


Equilateral Triangle 




A (0,0,0) 


(V3,0,l) 


(V3,0,-l) 




Isosceles - Apex Charged 




B (0,0,0) 


(0,0,2) 


(2,0,0) 




Line - Middle Charged 




C (0,0,0) 


(0,0,2) 


(0,0,4) 




Line - End Charged 




D (0,0,0) 


(0,0,2) 


(0,0,4) 




Isosceles - Leg Charged 




E (0,0,0) 


(0,0,2) 


(2,0,0) 




(b)Four sphere configurations. 


Configuration 


Letter Sphere A 


Sphere B 


Sphere C 


Sphere D 


Regular Tetrahedron 


A 


(y/2 V2 V2\ ( 
V 2 ' 2 ' 2 / v 


2 ' 2 ' 2 / v 


V2 V2 V2\ ( 
2 ' 2 ' 2 / v 


2 ' 2 ' 2 > 


Rhombus - Blunt Charged 


B 


(0,0,0) 


(0,0,2) 


(V3,0,l) 


(V3,0,3) 


Non-Reg Tetrahedron 


C 


(0,0,0) 


(2,0,0) 


(0,2,0) 


(0,0,2) 


Rhombus - Sharp Charged 


D 


(0,0,0) 


(0,0,2) 


(V3,0,l) 


(V3,0,3) 


Square 


E 


(0,0,0) 


(0,0,2) 


(2,0,0) 


(2,0,2) 



Appendix B: Explicit Matrix Equation 



This section of the matrix K (defined in equation 48) covers the projection of {Aq, A 1 1 , A®, A\} 
onto the monopole and dipole terms from spheres A, B, and the monopole from sphere 
C (given by vector {A^ A^.A^ A\, 5§, Sjf 1 , B\, C§}). For simplicity, we assume the 
spheres are all the same size [ta — — — 1)? an d we take the limit RaB) Rac, Rbc » 1- 
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This gives the first few rows and columns of IK as: 



K = 



47T 








0... 





47T 

3 





0... 








4-7T 

3 


0... 











47T 

3 ' * * 


47T 


2 v /fe-^AB 7 rSin[/3 A B] 


47rCos[/3 A B] 


2^/fe^AB 7rSin[/3 AB ] 


2y/l e ia AB 7rSin[/3 AB ] 


R 2 
n AB 

47rSin[/5 AB ] 2 


V3R 2 AB 
2V2e ia AB TrSinp^AB] 


/? 2 

n AB 

4e 2 ^AB 7rSin[/3 AB ] 2 


R AB 

4ttCos[/3 A b] 


3Ri B 

2V2e~ ia AB 7rSin[2/3 AB ] 


87rCos[/3 AB ] 2 


3i? AB 
2 v / 2e^AB 7rSin[2/5 AB ] 


Vsr 2 ab 

2y/l e— AB 7rSin[/3 AB ] 


3 4b 
4e -2za AB TrSin^AB] 2 


3R1b 

2v / 2e- iQ AB 7rSin[2^AB] 


3Rl B 
47rSin[/5 AB ] 2 


^AB 
47T 


3Rl B 

7rSin[/3 AC ] 


3-Ri B 

47TCOSI/3AC] 


3i? AB 
2^/fe^AC 7rSin[/3 AC ] 


RaC 


^AC 


V3R 2 AC 


R AC 



V 



/ 



As an example, the interaction between the m = dipole on sphere A and the m = dipole 
on sphere B is given by K 7 3 = — 87rCo ^ AB l since the seventh term in the vector M is £?? 
and the third term in M is A\. This term will also appear as K 3 j due to symmetry of the 
interaction. 



[1] When deriving the form for E BA , all odd order multipoles need to be multiplied by a — 1 
since the alignment of a positive multipole This adds a (— 1)^ factor to the odd ordered B 
multipoles in the derivation of E BA . 

[2] This makes the interactions much simpler by eliminating many of the azimuthal m ^ 
multipoles that occur in the general case interaction energy 

[3] Convention: E AB m is the interaction energy between a dipole on sphere A and a monopole on 
sphere B calculated with the E AB formula. E^ A q is the interaction energy between a monopole 
on sphere B and a quadrupole on sphere A calculated from the E BA formula 
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